wo 02/076585 PCT/US02/08706 



-1- 



ANALYSIS AND FRACTIONATION OF PARTICLES NEAR SURFACES 

RELATED APPLICATIONS 

5 This application is a continuation-in-part of U.S. Application Serial No. 

09/813,571 filed March 21, 2001, which is incorporated herein by reference in its 
entirety. 

FIELD OF THE INVENTION 

10 The present invention generally relates to the fields of analytical chemistry 

and materials characterization. Specially, the present invention relates to 
firactionation of heterogeneous mixtures of particles on the basis of characteristic 
frequency dependence of particle polarization and also to analysis of 
electrochemical properties of particles. 

BACKGROUND OF THE INVENTION 

LI Methods of Fractionation and Sorting 

Among analytical separation techniques, methods of flow cytometry have 
20 been used in sorting of particles (including cells) within a given heterogeneous 
population. These methods, however, require a dedicated instrument capable of 
selecting and physically separating from the population those particles which satisfy 
a given selection criterion. For example, the objective of cell sorting by flow 
2 J cytometry usually is the selection of those cells within a population displaying a 
characteristic fluorescence signal. 

L2 Methods of Particle Analysis 

Particle analysis represents a standard procedure of the analytical chemical 
repertoire that is used to determine physico-chemical and compositional properties 
[Hunter, "Introduction to Modem Colloid Science", Oxford University Press, 
Oxford, UK, 1993]. The extensive repertoire of techniques reflects the ubiquitous 
use of particles of a wide range of sizes, shapes, composition and chemical 
reactivity in many scientific and industrial applications ranging from medical 
35 diagnostics to cosmetics. Characterization of particles is useful in guiding and 
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o optimizing production as well as chemical modification particularly of surface 

properties that determine the stability of particle suspensions and the interaction of 
particles with molecules in the surrounding mediimi. 

Techniques of the standard repertoire may be grouped as follows: first, 

^ sedimentation and centrifiigation methods, electroacoustics, light scattering, 
hydrodynamic methods and dielectric spectroscopy; second, electrical pulse 
counting, flow cytometry and electrophoretic zeta potential measurements; third, 

spcuiaiiz.cu iiicuiuus liiuiuuui^ uicicuu^upiiuicsis {uE^r), iiiciuuiiig iiS curnDumuun 

with field-flow firactionation methodology. 
10 Methods in the first group apply to bulk suspensions and therefore deliver an 

average of the measured quantities for a large number of particles. While desirable 

in situations in which large industrial batches are to be characterized, bulk 

measurements generally do not lend themselves to miniaturization and are not well 
J ^ suited to the characterization of small numbers of particles, particularly when 

particle-to-particle variations are of interest or when multiple particle parameters are 

to be determined simultaneously. 

Methods in the second group use a "single-file" serial format of 

measurement. While these methods determine the properties of individual particles, 

20 

their implementation, in order to attain sufficiently high processing rates, requires 
high flow rates of a carrier gas or fluid and high-speed data recording and read-out 
electronics, generally rendering the equipment complex and expensive. This is so 
especially when multiple parameters are to be determined for each particle. For 
25 example, multi-color detection in flow cytometers can require the use of multiple 
lasers and multiple photomultiplier tubes. The determination of the zeta potential 
along with particle size in state-of-the-art equipment requires the measurement of 
electrophoretic mobility of particles traveling along a narrow capillary in response 
to a DC voltage applied along the channel as well as the application of light 
scattering with associated light source and read-out, as in the ZetaSizer (Malvem 
Instruments, Southborough, MA). Miniaturization, as in the case of flow cytometry 
[Fu, AX, Spence, C, Scherer, A., Arnold, F.H. and Quake, S.R. A Microfabricated 
Fluorescent- Activated Cell Sorter, Nature Biotechnology, Vol. 17, November 1999, 

35 
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o 1 1 09- 1 1 1 1 ] , retains the serial mode of processing and thus encounters the constraint 
of limited throughput. 

Methods in the third group include various dielectrophoretic techniques to 
characterize, classify and fractionate low conductivity particle suspensions [Becker, 

^ F.F., Gascoyne, P.R.C., Huang, Y. and Wang, X-B, Method and Apparatus for 
Fractionation Using Generalized Dielectrophoresis and Field Flow Fractionation, 
US Patent No. 5,888,370, March 30, 1999; Rousselet, J., Markx, G,H. and Pethig, 
R. , Separation of Erythrocytes and Latex Beads by Dielectrophoretic Levitation and 
Hyperlayer Field-Flow Fractionation, Colloids and Surfaces A: Physicochemical 

10 and Engineering Aspects, 140 (1998) 209-216; Pethig, R., Markx, G.H., Apparatus 
for Separating by Dielectrophoresis, US Patent No. 5,814,200, Sep 29^\ 1998; Pohl, 
H. A., Continuous Dielectrophoretic Cell Classification Method, US Patent No. 
4,326,934, April 27* 1982; Parton, A., Huang, Y., Wang, X-B., Pethig, R., 

J J MacGregor, A.R., and Pollard-Knight, D.V., Methods of Analysis/Separation, US 

Patent No. 5,653,859, Aug. 5, 1997; Crane, S., Dielectrophoretic Cell Stream Sorter, 
US Patent No. 5,489,506, Feb. 6^ 1996; Benecke, W., Wagner, B., Hagedom, R., 
Fuhr, G., Muller, T., Method of Continuously Separating Mixtures of Microscopic 
Dielectric Particles and Apparatus for carrying through this method, US Patent No. 

20 

5,454,472, Oct. 3^^ 1995]. All publications which are cited throughout the 
application are hereby incorporated by reference in their entirety. 

Coupled with a suitable imaging system, techniques within this group can 
provide single particle information. The particle classification is achieved by 

25 exposing the particle mixture to a non-uniform AC electric field which is generated 
by applying an AC voltage to multiple sets of planar patterned microelectrodes. 
Such a scheme in general requires a complex electrode signaling/addressing 
scheme. Also, a continuous mode of operation requires, in the simplest case, 
coupling to an external flow or the use of field flow fractionation. Simple 
dielectrophoretic setups do not generally allow for quantitative evaluation for the 
surface potential or surface conductivity of particles. While these quantities maybe 
determined from electrorotation spectra of single particles, this approach requires a 
complex and time-consuming set-up to confine a single particle at the origin of a 

3^ rotating electric field whose fi^equency is scanned while recording the rotational 
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o motion of the confined particle. Aside from the complexity of the equipment 

required to implement the requisite experimental configuration, the low throughput 
of the method presents a serious disadvantage when averages of multiple particles 
are desired. 

5 SUMMARY OF THE INVENTION 

The present invention provides a method and apparatus for the fractionation 
of a mixture of particles and for the analysis of particles, based in part on the 
technology designated "LEAPS", described in PCT International Application 
No.WO 97/40385 to Seul. The present invention uses LEAPS (referring to ''Light- 
ControUed Electrokinetic Assembly of Particles near Surfaces") to fractionate a 
heterogeneous mixture of particles on the basis of, e.g.: morphological properties 
including size and shape; surface composition relating to certain electrochemical 
properties, including surface potential and surface conductivity; and specific and 
non-specific adhesion to exposed surfaces provided within the apparatus. As used 
herein, the term "fractionation" refers to separation or sorting of pai:ticles, including 
separation of multi-component particle mixtures into their constituent 
subpopulations. Particles maybe separated on the basis of the characteristic 
20 frequency dependence of the particle polarization, said dependence reflecting 
physical or chemical properties of the particles. 

The present ihvention also provides for the use of LEAPS Xo determine 
various physical or chemical properties of particles, including morphological 
properties, surface electrochemical properties and surface chemical properties such 
as the particles adhesiveness for given surfaces of the apparatus or for other 
particles. 

The present invention provides a method and apparatus for fractionation of 
particle mixtures on the basis of the respective characteristic relaxation frequencies 

30 of the constituent particles or particle populations. In certain embodiments of such 
method, a plurality of particles are suspended at an interface between an electrolyte 
solution and a light-sensitive electrode. The particles comprise at least two types, 
each type having a distinguishable relaxation frequency. An electric field, having a 

2^ frequency that is less than or equal to the relaxation frequency of at least one of said 
particle types but greater than the relaxation frequencies of other particle types, is 
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o generated at the interface by application of an AC voltage. In addition, the interface 
is illuminated with a predetermined light pattern. The illumination in combination 
with the electric field produces fractionation of the particles having relaxation 
frequencies greater than or equal to the frequency of the electric field from other 
particles, e.g., by assembly of particles having a relaxation frequency greater than 
the frequency of the applied electric resulting in formation of a planar array of 
substantially one layer of particles in an area on the electrode designated by the 
pattern of illumination. In ccriaiii embodiments, the illumination affects the 
electrolyte-electrode interfacial impedance. 

10 In certain preferred embodiments, the electrode comprises a silicon electrode 

which is coated with a dielectric layer. In certain preferred embodiments, an 
additional electrode is provided such that the Hgjit-sensitive electrode and the 
additional electrode are substantially planar and parallel to one another and 
separated by a gap (e.g., in a sandwich configuration), with the electrolyte solution 
containing the particles being located in the gap. 

In certain other embodiments of the present invention, patterned electrode is 
used instead of or in combination with the illumination to fractionate a mixture of 
particles. In one such example, a first electrode positioned in the first plane and a 
second electrode positioned in the second plane different from the first plane is 
provided with a gap between the electrodes , an arrangement herein referred to as a 
sandwich configuration. A mixture of particles is suspended in an electrolyte 
solution and located in the gap. The second electrode comprises a planar electrode 

25 having a surface and an interior, either of which may be patterned to modify the 

spatial distribution of the interfacial (interface between the electrolyte solution and 
the second electrode) electric field. When an electric field is generated between the 
first and the second electrode by applying an AC voltage between the two, the 
electric field in combination with the patterning of the electrode produces 

30 

fractionation of the particles according to relaxation frequency; for example, 
particles having a relaxation frequency greater than the frequency of the applied 
electric field move toward designated regions of the second electrode surface as 
provided by illumination or patterning and therein form a planar array of 
35 substantially one layer. In preferred embodiments, the second electrode may be 
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pattemed by spatially modulated oxide growth, surface chemical patterning or 
surface profiling (which is described in U.S. Serial No. 09/171,550, filed October 
26, 1998, now U.S. Patent 6,251,691 issued June 26, 2001); preferably, the 
patterning affects the impedance of the electrolyte solution and second electrode 
interface. In preferred embodiments, the second electrode comprises a silicon 
electrode which is coated with a dielectric layer. 

In certain embodiments, the difference in relaxation fi*equency of particles 
within a heterogeneous mixture corresponds to a difference in particle size or 
surface composition, thus permitting fractionation of the mixture into homogeneous 
subpopulations of particles according to size or surface composition, respectively. 
The process of firactionation may also be monitored using a video detector or 
camera. 

Accordingly, the present invention provides a method and apparatus for 
fi-actionation of particle mixtures based on the characteristic relaxation frequencies 
of said mixture's constituent particles or particle subpopulations. The present 
invention provides for the determination of the frequency dependence of the 
electrohydrodynamic forces that govern particle transport into designated low- 
impedance areas of the electrode and mediate formation of particle arrays in those 
designated areas. As described herein, a particle relaxation frequency, cor, is a 
property that is associated with the particle-solution interface and may be 
determined from the condition that array assembly requires co<cor,: particles not 
meeting this criterion will be excluded from low-impedance areas. 

The present invention also provides for a method and apparatus for the 
analysis of particles in order to determine various physical or chemical properties of 
said particles, including size, morphology, surface potential and surface 
conductivity. For example, the relaxation frequency and/or maximal velocity of 
transport of the particles are used for particle characterization. 

In certain embodiments of the present invention, a plurality of particles 
suspended in an interface between an electrolyte solution and a hght-sensitive or a 
patterned electrode are provided, as described above in connection with particle 
fractionation. When an electric field is generated at the interface by the application 
of an AC voltage, the firequency of the electric field is adjusted to produce formation 
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o of a planar array of substantially one layer of particles in an area designated by the 
pattern of illumination or electrode patterning. From the analysis of particle 
transport in the course of array assembly, the relaxation frequency and the maximal 
velocity of transport of the particles may be determined. AUematively, the particle 
relaxation frequency also may be determined by analyzing the response of the array 
configuration to changes in the frequency of the applied electric field. 

The present invention also provides for the determination of the frequency- 

denendent and vn\f^^xffi^-^^p1t^p>J\f^ev^i maYimsil vplnriH/ V . aftQin^H hi/ noH-iVlAc 
crossing impedance gradients in the course of array assembly. As described herein, 

*0 Vmax may be determmed by image analysis and particle tracking. In certain 
embodiments, the combination of cor and Vmax provides for the simultaneous 
determination of the particles' surface ("zeta") potential, (|)s, and surface 
conductivity, c, e.g., based on numerical analysis of phenomenological equations of 

15 particle motion described herein. 

Once these values are found, they may be used to characterize the particles, 
e.g., determine the zeta potential of the particles and/or the mobility of ions or 
molecules within the electrolyte solution in proximity to the particle surface. In 
preferred embodiments, the relaxation frequency and the maximal velocity are used 
together in particle characterization. In certain preferred embodiments, the present 
invention permits the simultaneous determination of multiple particle properties, 
such as zeta potential and ionic or molecular mobility. In other embodiments, the 
present invention permits the determination of one particle property of interest 

25 independently from the determination of others. For example, in certain 

embodiments, the ionic or molecular mobility associated with the particle-solution 
interface may be determined from the relaxation frequency alone. 

The term "particles," as used herein include colloidal beads, eukaryotic and 

3Q procaryotic cells, micelles, vesicles, and emulsion droplets. In preferred 

embodiments, the particles comprise colloidal beads, or eukaryotic or procaryotic 
cells. 

In contrast to prior art methods of fractionation of particle mixtures such as 
field flow fractionation, the methods of the present invention provide a parallel 
format based on the manipulation of a multiplicity of individually imaged particles 
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in a simple sandwich geometry which provides for longitudinal and transverse 
electrohydrodynamic forces to mediate fractionation on the basis of the frequency- 
dependent polarization of the particles of interest. Specifically, the methods of the 
present invention provide for the fractionation of one subpopulation of particles at a 
time from the remainder. 

In contrast to prior art methods of sorting such as fluorescence-activated 
particle and cell sorting, the methods of the present invention provide for a parallel 
process accommodating a multiplicity of particles thereby enhancing processing 
speed. The methods of present invention also are well suited for miniaturization. 

Li contrast to prior art methods of particle analysis which are applied to bulk 
particle suspensions and yield ensemble averages of the particle properties of 
interest, the methods and apparatus of the present invention provide for the 
miniaturization of the processes of analysis. In contrast to prior art serial methods of 
analysis of particle mobiUty, the methods of present invention invoke a collective 
phenomenon of particle array assembly as the basis for a parallel mode of 
simultaneous analysis of multiple particle properties accessible either to imaging 
and image analysis or to light scattering. The methods of the present invention 
provide for the simultaneous, yet independent determination of morphological 
parameters including size, orientation and polydispersity, as well as for 
electrochemical parameters including surface potential and surface conductance. 

In contrast to these prior art methods, the implementation of the fractionation 
methods of the present invention are readily suitable for miniaturization and 
integration with additional preparative and analytical procedures in an integrated 
analytical environment. 

BRIEF DESCRIPTION OF DRAWINGS 
Fig. 1 is an illustration of the response of a polarizable double layer around a 
spherical particle to a time-varying electric field. The progression of polarization 
states occurring in a single cycle is shown in panels (a) - (e). The shading indicates 
the distribution of the predominant type of ion in the diffuse layer. 

Fig. 2 is an illustration of the variation of the polarization factor of a particle 
with frequency, as evaluated numerically on the basis of the theoretical model 
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described herein and assuming the indicated values of the particle surface potential, 
a particle diameter of 2.8 \im and a solution of ionic strength 0. 1 mM. 

Fig. 3(A) is an illustration of the dependence of the relaxation frequency on 
surface potential, as evaluated numerically on the basis of the theoretical model 
^ described herein for a 2.8 ^lm particle suspended in a solution of ionic strength 0. 1 
mM. The legend refers to the ratio of counterion mobility in the double layer to the 
bulk ion mobility. (B) is a plot of the relation between surface potential and 
mobility for the specific relaxation frequency of 1000 Hz, 

Fig. 4 is an illustration of flow fields characteristic of a sandwich geometry 
10 compatible with LEAPS. 

Fig. 5 is a flow chart summarizing an analytical procedure for the 
simultaneous determination of particle surface potential and surface conductance 
from experimentally determined values of particle velocity during array assembly 
12 and particle relaxation frequency. Also indicated are the computer programs 
("Codes") generated to perform the calculations. 

Fig. 6 is an illustration of particle movement and of the relative position of 
each particle relative to a position defined as x=0. The result of such particle 
movement is array assembly and each particular particle position represents a stage 
in the assembly process. The particles assemble in thin oxide region of the 
silicon/silicon substrate and the demarcation between thin oxide and thick oxide 
defines the position x=0. 

Fig. 7(A) is an illustration of the particle velocity at the interface as a 
25 function of surface potential, as evaluated numerically on the basis of the theoretical 
model described herein in for several ratios of the counterion mobility to the bulk 
ion mobility, as indicated in the legend. (B) is a plot of the relation between surface 
potential and mobility for the specific bead velocity of 0,8 mm/s. 

Fig. 8 is an illustration of the process particle fractionation by selection of 
frequency, co (a)R^< o) < cdr^) from an initially random binary mixture of particles 
having two different relaxation frequencies (cdr^ and o)r^, cor^ > cor^) using a 
substrate surface featuring a half plane of lower impedance adjacent to one of higher 
impedance. 

35 
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o Fig. 9 is an illustration of the determination of the velocity of particles 

experiencing acceleration as they cross a boundary, located at position "0", 

separating regions of high and low electrode impedance. 

Fig. 10 is an illustration of experimentally determined dependence of 
^ particle relaxation frequency on particle diameter for three different types of 

particles. Experimental details are indicated in the inset to the graphical display. 

The results for three different saU concentrations are shown for two ionic diffusion 

coefficients. Result.^ show that th<» rAlnvsitinn fri»ni 

■ ' " " — — ■ Au &«v/« 0&^ui.iwcui^ijr c&xxv«v/iV/U 

by a change in salt concentration. 

Fig* 1 1 is an illustration of the process of fractionation of an initially random 
particle mixture composed of larger and smaller particles on the basis of particle 
size using a substrate surface featuring a strip of lower impedance (bright area) 
surroimded by areas of larger impedance (dark areas). 
^5 Fig- 12 is an illustration of the two successive stages in the process of 

fractionation of an initially random particle mixture composed of particles of 
identical size which display either the protein streptavidin (shown in red) or 
displaying a layer of short oligonucleotides attached to the streptavidin-coated 
particles by way of a biotin functionality (shown in turquoise). 

Fig. 13 is an illustration of the collection and assembly of HL-60 Human 
Myeloid Cells (-15 ^im diameter) in characteristic configurations on a patterned 
silicon chip. 

Fig. 14 is an illustration of the collection and assembly of E. coli and staph 
25 bacterial cells (-1 \im diameter) in characteristic bands on a patterned sihcon chip. 

Fig. 15 is an illustration of serial fractionation of a multi-component particle 
mixture (with each particle type possessing a unique relaxation frequency(o)R )using 
patterned electrode and sequence of illumination patterns. 

30 



20 



35 



wo 02/076585 



PCT/US02/08706 



10 



15 



- 11- 

DETAILED DESCRIPTION OF THE INVENTION 

I. ELECTROSTATIC THEORY OF DIELECTRIC MEDIA 

So that the present invention may be better understood, certain results from 
the electrostatic theory of dielectric media are herein summarized, making general 
reference to Landau, Lifshitz & Pitaevskii, "Electrodynamics of Continuous 
Media", Butterworth-Heinemann, Oxford, UK, 1998. 

I-l PoDderomotive Forces 
The free energy of a dielectric medium of dielectric constant tu filling a gap 
of size H between two parallel capacitor plates of area A carrying a charge Q is 
given by 3o = -VzCq^^, where the capacitance, Co = em A/H, describes the 
proportionality between applied potential and stored charge, Q= Co(t). Under the 
assumption that the electrostatic potential of the plates remains unchanged, the 
introduction of a sphere of radius a « H and dielectric constant into the gap 
changes the free energy in accordance with the following expression: 



20 



25 



30 



(1) 



Here, the factor * = (Ss - emVCes +2eM) determines the potential of an 
uncharged dielectric sphere of radius R and dielectric constant embedded in a 
dielectric medium of dielectric constant em and subjected to a homogeneous electric 
field Eo. 



^= -£^>rcos^ 








r>a 









(2) 



In accordance with Eq.(l), as long as eM< Cs, insertion of the sphere into the 
dielectric medium in fact lowers the free energy. Consequently, a sphere placed just 
outside the edges of the plates will be drawn into the space between the plates by a 
force V3. From Eq. (1), with £o = <t> 



F = Ine^e^kVE' 



(3) 
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This ponderomotive force thus arises from the interaction of the polarization 
of the sphere with spatial non-uniformities in the electric field. Eq.(3) shows that the 
force is independent of the sign of the electric field. That is, in a non-uniform 
electric field, a polarizable particle will move either into or out of a region of higher 
field, as determined solely by the relative magnitude of the dielectric constants of 
particle and suspending medium. In general, as discussed herein, (he polarization 
factor k will deoend on fi-eouencv! k = hfff\\ 

1.2 Interfacial Polarization in Electrolyte Solutions 

The polarization of a fluid medium, herein referred to as an electrolyte 
solution, typically manifests itself in the form of gradients in the spatial distribution 
or orientation of the constituent molecules. Thus, electrolyte solutions are polarized 
near any interface with a bounding (boundary) surface or interface. When such an 
interfacial region is subjected to an external electric field of varying fi*equency, the 
15 response of the medium is determined by the processes mediating redistribution or 
reorientation of the constituent molecules and by characteristic molecular mobilities 
associated with these processes. 

Spatial Distribution of Free Ions: The Double Layer - The response of an 
2Q electrolyte containing mobile charges to the presence of a charged bounding surface 
is the formation of a layer of excess charge adjacent to the surface so as to screen 
the field generated by the charges on that surface. Thus, the polarization of the 
interface between a charged particle and an electrolyte containing free ions 
manifests itself in the form of a redistribution of ions so as to form screening shell 
of excess counterions. This "double layer" is characterized by a surface potential, ({>s, 
which corresponds to an excess of free charge carriers adjacent to the sphere and by 
a characteristic screening length, k\ which is directly related to the ionic strength, 
!/2ZnioZi^ of the medium; for a 1 :1 electrolyte, k * = (eMeokBT/2z^e^no)^ where no is 
30 the molar concentration of the electrolyte and z is the charge on each ion. The local 
excess concentration of ions can affect binding reactions involving molecules 
anchored to the sphere surface. A uniform electric field induces the redistribution of 
ions adjacent to and within the "double layer", leading to a phenomenon sometimes 
2 J referred to as "concentration polarization" (Fig 1). The redistribution of the 
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o molecules forming the screening layer is determined by translational molecular 
mobilities. 

Preferred Orientation of Dipolar Molecules - The response of an electrolyte 
containing dipolar molecules to the presence of a charged bounding surface is the 
formation of a layer in which those molecules assume a preferred orientation. As 
with the previous case, the interfacial region of preferred molecular orientation is 
characterized by a surface excess polarization and a characteristic decay length 

which is readilv HeterminftH frnm p ctfltictiml tViArrt-\/ Th** Ir^nal 

J — — w>ww> J , A SAW W^WWkJO Vm AWUMAWIVIA 

can affect binding reactions involving molecules anchored to the sphere surface. 
JO The redistribution of the molecules forming the screening layer is determined by 
orientational molecular mobilities. In any given electrolyte, excess concentration 
and orientation (as well as higher moments of the molecular charge distribution) 
may contribute to the formation of a screening layer. 

15 1.3 Frequency Dependence of Field-Induced Polarization 

An externally applied electric field alters the interfacial polarization of 
electrolytes by inducing the redistribution or reorientation of the constituent 
molecules. Thus, an applied AC electric field induces a time-dependent interfacial 
2Q polarization. As with the recharging of a parallel plate capacitor, the time-dependent 
interfacial polarization corresponds to a transient charging current,/ ~ ft When 
the time-dependent interfacial polarization is mediated by the spatial redistribution 
of free ions, the charging current,/, will be composed of a longitudinal and a 
transverse component, respectively reflecting the redistribution of ions within the 

25 

double layer and the exchange of charge between the double layer and the 
surrounding medium. Set forth below are certain results for this case; it is noted that 
corresponding results also may be cited for the time-dependent reorientation of 
dipolar molecules. 

30 The longitudinal component reflects the redistribution of primarily of 

counterions within the double layer in response to the applied field on a 
characteristic time scale T|| - (a + 1/k)^/Z) - I/cor, where a, 1/k and D respectively 
. denote: the particle radius, the Debye screening length and a characteristic 

2^ dififusivity of the majority carrier. The characteristic time scale, x = X[|, is determined 
primarily by the counterion mobility, u - ze/f^ where, n and / respectively denote the 
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local counterion concentration and drag coefficient. Ionic mobilities in the vicinity 
of the particle surface will reflect the details of the local chemical environment and 
generally will differ significantly from the corresponding value in the bulk medium. 
By invoking the Einstein relation \lf= Dlk^T, the relaxation frequency, (wr, may be 
expressed in terms of the counterion mobiUty, u: 




(4) 

The longitudinal counterion current is balanced by a transverse current 
reflecting the radial transport of co- and counterions across the double layer on a 
time scale xi \I\^D 1/coi. If electroneutrality is to be strictly maintained in the 
electrolyte adjacent to the double layer, then the exchange of counterions between 
double layer and surrounding electrolyte via radial transport must be accompanied 
by the radial transport of an equal number of coions. The result is a change m the 
local ionic strength in the electrolyte adjacent to the double layer. That is, gradients 
in electrolyte concentration appear on the scale of the particle diameter, and these 
decay away on a characteristic time scale x - 1/cor. As the frequency increases, the 
20 redistribution of ions within the double layer eventually will no longer suffice to 
attain the fully polarized state; relaxation occurs when ca^x = 1 . For a typical 
particle radius of 1 ^m, and a typical value cOj^ - 10^ /s, the relaxation frequency, 
CO-- coi/(l + Kfl)^ is approximately 1 kHz. 

25 

Phenomenological Theory for Polarization Factor k(6}) - Phenomenological 
models have been developed to relate the field-induced, frequency-dependent 
polarization, k = /:(co) to the particle's surface potential, (|)s, and to the surface 
conductivity, a. The former, for given bulk ionic strength, determines the total 

30 charge density within the double layer; the latter determines the rate at which the 
double layer is able to readjust in response to a time-varying applied field. 
Specifically, a linearized theory, valid for small surface potentials, is available to 
treat the contribution of transverse ionic fluxes to the transient ionic charging 

2 J currents and the charge balance within the double layer to produce an expression for 
the frequency dependence of A: = /:(a)) [Dukhin and Shilov, "Dielectric Phenomena 
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and the Double Layer in Disperse Systems and Poiyelectrolytes", Keter Publishing, 
Jerusalem, 1974]. A recent formulation of the linearized theory of the electrokinetic 
transport phenomena underlying concentration polarization, developed for the 
analysis of the low-frequency dielectric response of colloidal suspensions, yields a 
useful analytical form for the frequency dependence of A<o) [Minor et al, J. Colloid 
& Interface Science 206, 397-406 (1998)]. 

Assuming low values of the external field, the potential and the ion 

UAtJU jLri4.t^X\/Aiyi9^ llAAVUI.t£/VU CUVfUIlU. LAAV V\|lUXlLri lUlll WllUlilVril VA.i3VlXi^ ill 1111/ 

absence of an external field. For the region within the suspending medium far from 
the double layer, ion balances (for the i-th species) obey the Nemst-Planck equation: 



Here, the terms preceded by *5' denote the respective first order deviations 
fi*om the corresponding equilibrium values, and f is an effective fiiction factor 
which governs the ion mobility. Solving Eq. (5) along with the linearized Poisson 
equation, expressions for 5ni ^ and fare found as follows: 



Here, Xa = (1+z) (fflx)^; the time scale, x, is set by the characteristic dififiision time 
for transport aroimd the perimeter of the particle: r= a^/D -f*c?/2k^T 

The polarization parameter, K determines the increase in the effective 
dielectric constant of the medium due to the presence of particles and associated 
double layers. To express k = /:(co) in terms of system properties of interest, 
boundary conditions are applied to Eq.(6) in the form of a flux balance: tangential 
ion fluxes within the double layer are balanced by radial ion fluxes into and out of 
the double layer so as to ensure that electroneutrality is maintained in the external 
solution at all times. Restricting attention to an electrolyte containing only a single 





wo 02/076585 



PCT/US02/08706 



10 



15 



20 



-16- 

coion (superscript "1") and a single coimterion species (superscript "2") and 
ignoring the coion flux within the double layer as negligible compared to the 
counterion flux, the flux balance at the outer edge ( r - a"*) of the double layer is 
expressed in the form: 



+ Du2aV] 



= 0 



(7) 



where V, denotes a tangential gradient operator. Radial and tangential counterion 
fluxes are related by the constant, Du2, which represents the dimensionless ratio of 
(the counterion-mediated) surface conductivity to (the counterion-mediated) bulk 
conductivity, Du2 = a/aaoo, a denoting the particle diameter. Du2 may be expressed 
in terms of the surface potential, ^s, and a dimensionless counterion mobility, /W2, as 
follows: 



, where wj, = 



2e,kTf,' 



(8) 



Application of the boundary conditions yields the following expression for 
the polarization parameter, k{(a): 



1 3 

Kw = + f 5, 

2 2 2 + Duj[2-y(l-h)] 



where 



Y^l- 



1 + Xa 



COT 



+ ICOT 



and 



Fig. 2 shows representative plots of A(co) for different values of surface 
potential and for a given ratio of surface ion mobility, w, and bulk ion mobility, 
The point of inflection in each of these plots identifies a peak in the dispersion. 
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o which in turn identifies the corresponding frequency as the relaxation frequency, 
cor. The plots in Fig.3 show that for given value of cor generally will exhibit a 
dependence on surface potential. That is, cdr may be understood as a function of the 
ionic conductivity, a eun ze^n/f^ a quantity which in turn depends on both ionic 

J mobility as well as local ion density, the latter being controlled by the surface 
potential. 

However, Fig, 3 also shows that, for sufficiently large values of the surface 
potential, <{)$, the dependence of cdr on surface potential vanishes, yielding instead a 
linear relationship, cor~ u/u^, in accordance with the scaling relation, Eq. 4. This 

10 

relation, which formed the basis for the analysis of the data in Fig. 10, thus provides 
the basis for a direct experimental determination of ionic or molecular mobility from 
a determination of the relaxation frequency. 

II. ELECTROKINETIC PARTICLE TRANSPORT 
AND PARTICLE ARRAY ASSEMBLY 

Light-controlled Electrokinetic Assembly of Particles near Surfaces 

(LEAPS), as disclosed in W097/40385, enables the formation of planar arrays of 

particles (including on-cue formation) in designated areas of a substrate under real- 

20 time optical control. LEAPS combines the use of electric fields, interfacial 
patterning and/or light to introduce spatial and temporal modulations of the 
impedance of an Electrolyte-Insulator-Semiconductor (EIS) interface. LEAPS 
invokes electrohydrodynamic forces which arise in accordance with the lateral 
impedance gradients to control the coupled transport of fluid and suspended 
particles as well as the array assembly process. 

The present invention describes the role of polarization-induced forces that 
contribute to particle transport and array assembly under conditions characteristic of 
LEAPS. These models guide the selection of conditions for the fractionation of 

30 particle mixtures in response to combinations of time-varying electric fields. These 
models also provide methods for the analysis of coupled fluid-particle transport in 
order to determine characteristic particle properties of interest. 



35 
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o IIL FLUID FLOW, PARTICLE COLLECTION AND 
ARRAY ASSEMBLY 

Light-controlled Electrokinetic Assembly of Particles near Surfaces 
(LEAPS) is described in PCT International Application No. PCT/US97/08159 
(published as WO 97/40385) and U.S. Serial No. 09/171,550, filed October 26, 
5 1998 (issued as U.S. Patent 6,251,691 on June 26, 2001) which are incorporated by 
reference in their entirety. Also incorporated herein by reference in their entirety are 
U.S. Serial No. 09/397,793, filed September 17, 1999 (pubHshed in WO 01/20593) 
(describing system and method for programmable illumination pattern generation), 

jQ as well as U.S. Serial No. 09/320,274, filed May 28, 1999, and PCT International 
Application No. PCT/USOO/14957 (published as WO 00/73777) (relating to array 
cytometry). LEAPS involves the interplay between several electrokinetic forces to 
induce fluid flow, particle transport and array assembly. For example, in a sandwich 
geometry, an AC electric field applied across a fluid gap generates rectified 
electroosmotic flow in the vicinity of lateral gradients of the Electrolyte-Insulator- 
Semiconductor (EIS) impedance that are set up by interfacial patterning or spatially 
varying illumination (see Fig. 4). Generally, these flows are observed under the 
condition (o < ©i, the latter characteristic fi-equency being determined by molecular 

20 processes governing the response of the polarization of the EIS interface. Under the 
action of these flow fields, particles are "collected" over macroscopic distances and 
are swept from regions of higher toward regions of lower interfacial impedance by 
viscous drag. 

22 Polarizable particles, while being swept along by the electrokinetic flows, 

simultaneously respond to the applied electric field; specifically, if particles are 
charged, the applied field induces local charging currents which in tum create a 
substantial concentration polarization in accordance with the mechanism previously 
discussed herein. This field-induced polarization determines the effective dielectric 

30 

constant of polarizable particles and so dominates the ponderomotive forces exerted 
on such particles in regions of rapid spatial variations in electric fields such as those 
produced near steps or other steep gradients in interfacial impedance. 

Accordingly, particles approaching a step in EIS impedance come under the 
35 influence of the gradient V{E'E) and accelerate across the step into the low 
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o impedance region where they are slowed, e.g., by the drag force exerted by the 
fluid. The particles enter into regions of low impedance and form airays when 
CO ^ Or ("aggregation"); in contrast, when co > cdr, particles accumulate on the high 
impedance side in the vicinity of the boundary. Illumination also can induce 

^ gradients in impedance and, provided that © ^ ©r, array formation in illuminated 
areas of the electrode. 

In certain embodiments of the present invention, the frequency dependence 

fit A rtm-wn-mr -fV^wM « A-n «%«>Ar«rkO«« MAwMvfn fl>r« Aii4 1 ^ Af a-r-rt^^v«rtf 't AM <%f fit A 

characteristic frequency, or. For a homogeneous population of particles, the 
^® determination of the relaxation frequency is determined by inducing a transition in 
the array configuration: particle arrays assembled under the condition © < cor 
assume a close-packed configuration; when the frequency is then tuned to co > cor, 
arrays undergo a transition to an expanded configuration [Chau et al, Proc, SPIE, 
15 3877, 36-42] . Accordingly, the frequency of the applied AC electric field is first set 
to a sufficiently low value so as to enable the assembly of an array within a 
designated area and is then set to progressively higher values to trigger the transition 
to an expanded array configuration. 
2Q For a heterogeneous mixture of particles disposed on a pattemed electrode 

surface, the entire set of relaxation frequencies characterizing the types of given 
particles maybe determined as follows. The frequency of the applied AC electric 
field is first set to a sufficiently high value so as to prevent transport of any particles 
into designated areas of low impedance, and is then set to progressively lower 

25 

values to permit array assembly for those particles satisfying the criterion © ^ ©r. 

Particle Size - Using this transition in array configuration to determine ©r, 
the relation between ©r and particle size may be determined for different electrolyte 
concentrations and for particles displaying different surface properties (see Example 
30 2). 

Surface-Chemical Composition - For given size, ©r also varies with particle 
surface chemistry. For example, when proteins are coupled to bead surfaces, ©r 
decreases; when oligonucleotides are coupled to protein-coated beads, ©r, increases. 
35 In general, both the local density of charge carriers, «, as well as their ionic 
mobility, w, may be affected. 
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o Non-aqueous Electrolytes - Rectified field-induced fluid flow as well as 

particle transport and array assembly, the latter accompanied by an acceleration 
across impedance steps, also occur in non-aqueous suspensions. The relaxation 
frequency of particles assembled in DMSO and acetonitrile is shifted to lower 

^ values when compared to aqueous suspensions. 

IV. FRACTIONATION OF HETEROGENEOUS PARTICLE 
MIXTURES 

The characteristic dependence of oor on panicle size and surface chemical 
composition, manifesting itself in its effect on particle transport and array assembly 
under conditions characteristic of LEAPS as described herein, provides the basis for 
an on substrate (e.g., on-chip) method to fractionate complex heterogeneous 
mixtures of particles including suspensions of cells. The underlying principle of 
operation of the method of the present invention is to form designated areas of low 
15 impedance of an EIS interface, either by interfacial patterning or by illumination as 
provided in LEAPS, and to select a frequency, or sequentially select multiple 
frequencies, of the applied voltage, under conditions compatible with the 
aggregation of at least one constituent subset among a given, heterogeneous set of 
particles. 

20 

To fractionate a mixture composed of two different types of beads that differ 
in their respective relaxation frequencies, the frequency of the applied voltage is set 
to an intermediate value to selectively assemble arrays of beads of the higher 
relaxation frequency. In accordance with the dependence of relaxation frequency on 
particle size, as disclosed herein, the condition (Or'^< co < cor^ ensures that only the 
smaller particles (characterized by a frequency (Or) enter the region of lower 
impedance, while larger particles (characterized by a frequency cor^) align on the 
high impedance side of the feature (see Example 4). Similarly, the dependence of 

30 the relaxation frequency on surface chemical composition can be used to fractionate 
particles (see Example 5). 

The binding of ligands to the surface of fimctionalized beads that produces a 
change in surface chemical composition can be detected on the basis of the 

2^ corresponding shift in relaxation frequency. In its simplest format, the process is 
illustrated by the selective aggregation in a low impedance area of the electrode of 
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only one of an originally randomly mixed set of two types of particles, with 
respective relaxation frequencies ©r^ < cor^, by selection of the operating frequency, 
©R^ < (Or< ©r^, as illustrated in Fig. 8. 

This fractionation process of the present invention may be generalized to 
multi-component mixtures by progressively selecting frequencies to isolate and 
assemble one constituent subset of beads. The method of the present invention also 
applies to cells whose differential adhesion properties provide an additional 
discriminating feature for LEAPS-mediated fractionation. By the method of the 
present invention, different types of bead and cells are fractionated and sorted 
without the need for complex and expensive equipment such as flow cytometers. 

As with the method disclosed herein, polarization-dependent forces also play 
a role in dielectrophoretic methods for fractionation of heterogeneous cell 
populations. However, the methods of present invention provide several distinct 
advantages, including: 

- the combination of electrokinetic flow and polarization- 
dependent particle responses in a collective process of particle array assembly, 
thereby providing for the integration of particle transport with fractionation in 
designated areas of an electrode; 

- fractionation on the basis of surface chemical properties, thereby 
providing for the detection of changes in the surface composition of particles or 
cells in the course of binding assays; 

- sequential multi-step fractionation within substantially the same 
area of a substrate under optical monitoring; 

- fractionation by multiple parameters including size, surface 
composition and selective adhesion; 

- multi-dimensional fractionation, e.g., 

frequency dispersion of array assembly, followed by selective adhesion; 

- parallel, multi-contrast mode of detection via optical imaging; 

- ease of use due to simple handling of electric contacts and; 
flexible electrode patterning by standard methods of semiconductor processing; 

- compatibility with dielectrophoretic separation processes; 
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- optical programmability as provided by Seul and Chau [Chau et 
al» Proc. SPIE, 3877, 36-42]. 

Several related distinctions may be drawn between the method of the present 
invention and dielectrophoretic methods of particle separation including: 
^ Underlying Process - The concentration polarization process invoked by the 

present invention, detectable by dielectric spectroscopy of particle suspensions 
[Minor et al, J. Colloid & Interface Science 206.397-406 (1998)], probes molecular 
configurations and transport phenomena associated with particle surfaces. For 
particles within the range of sizes of typical interest to the applications considered 
^0 herein (0.1 ^m - 10 \xm), this process is associated with characteristic relaxation 
frequencies in the range of a few hundred Hz to a few tens of kHz, a range 
comprising AC field-driven electrokinetic flow phenomena associated with the 
electrode surface [Yeh et al. Nature 386, 57-59 (1997)]. 
15 In contrast, representative implementations of dielectrophoretic methods for 

fractionation of particle mixtures developed by Pethig's group [G. H. Markx, M, S. 
Talary, R. Pethig, J. Biotechn 32, 29-37 (1994)] and by Gascoyne's group [Y. 
Huang, X-B. Wang, F. F. Becker & P. R. C. Gascoyne, Biophys. J. 73, 1118-1 129 
(1997)] use the fi-equency-dependent transition from a positive to a negative 
dielectrophoretic response of cells at high frequency (typically exceeding lOOkHz) to 
distinguish between cell types; this transition, detectable by the response of 
individual cells to rotating electric fields in electrorotation measurements, reflects 
the relative contributions of the dielectric permittivity and membrane conductance 
25 to the overall dielectric response of cells [Jones, "The Electromechanics of 

particles", Cambridge University Press, Cambridge, UK 1995]. In contrast to the 
collection of particles over macroscopic distances afforded by electric-field-induced 
fluid flow in accordance with the present invention, the prior art relating to 
2Q dielectrophoretic methods of particle fi-actionation does not invoke electric field- 
induced fluid flow in the fractionation process. 

Electric Field Configuration - As with LEAPS, the configuration of electric 
fields according to the method of the present invention is produced by applying an 
AC voltage normal to the EIS interface, for example by adopting a sandwich 
configuration in which the voltage is appUed between a pair of parallel planar 
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o electrodes confining the electrolyte solution. That is, the apphed field is nominally 
normal to the electrode surfaces. In contrast to the methods of the present invention, 
voltages in these dielectrophoretic methods of separation are not applied across the 
fluid gap, but are instead applied between pairs of electrodes deposited on the same 
^ surface. Because the forces of interest in these methods of the prior art are of very 
short range, arising primarily firom the spatial inhomogeneities of the "fiinge" fields 
at the electrode edges, designs typically employ sets of closely spaced metal 
electrodes of opposite polarities, 

Pethig and colleagues use interdigitated, "castellated" microelectrode 

10 designs and scan frequency to separate cells by positive and negative 

dielectrophoresis; for example, a range of fi-equency can often be found m which 
viable cells experience positive dielectrophoresis while non-viable cells experience 
negative dielectrophoresis. Differential forces of retention strongly confine cells 

J J under positive dielectrophoresis, permitting cells experiencing negative 

dielectrophoresis to be swept away by an externally provided lateral fluid flow. 
Gascoyne and colleagues have devised a dielectrophoretic levitation method in 
which traveling wave dielectrophoresis is combined with externally driven laminar 
flow along the electrode surface in a field-flow fractionation process in which the 

20 

differential dielectrophoretic lift of different types of cells exposes them to different 
flow velocities of a parabolic flow profile maintained in the separation chamber. 
Entrained by the flow, cells are detected, one at a time,. as they traverse a 
downstream cell counter. 

Simultaneous Transport and Assembly - The LEAPS-enabled fractionation 
method of the present invention combines electrokinetic flows with concentration 
polarization in a collective process to integrate flow-mediated particle transport with 
the differential response to ponderomotive forces to bring about both the long-range 
transport of particles ("collection") as well as fi-actionation in the process of array 
assembly. The fi-equency ranges typically employed in dielectrophoretic analysis 
exceed the characteristic relaxation ftequencies associated with electric field- 
induced electrokinetic flows, a fact which has profound consequences for the design 
of electrode configurations permitting the short-ranged dielectrophoretic forces to 
be utilized for particle collection. Voltage switching strategies (as for example in 
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traveling wave dielectrophoresis), enabling serial transport and fractionation of 
particles, however, suffer from the disadvantage that it requires complicated 
electrode configurations and design. 

Optical Progranmiability - The method of the present invention is unique in 
providing for real-time control of fractionation processes by means of light-induced 
interfacial impedance modulations afforded by the programmable generation of 
illumination patterns. 

V. PARTICLE ANALYSIS 

The present invention includes a method of analyzing particles based on the 
analysis of particle transport in the course of particle array assembly in accordance 
with LEAPS. To that end, the equation of motion of a spherical particle including 
viscous drag and ponderomotive forces may be written in the form: 



|toV^= -6;r/ia(v,- Vf)+ ^(<y)VE^ 



(9) 

where Vb and Vf respectively denote the velocity of the particle and the fluid medium 
and K((o) = 27ia^SMA, where k is defined in connection with Eqs (1) and (2). The 
ponderomotive force will vanish as A: -» 0 or Ss Bm- Reference to Fig. 2 shows 
that this condition corresponds to relaxation: as the frequency of the external electric 
field approaches cor, the fiilly polarizable state is no longer attainable. In general, 
Vmax will thus depend on surface potential, ion mobility and frequency. 

This equation may be used to analyze the motion of particles in the course of 
particle collection and array assembly by implementing the following steps: 

1 . For a given BIS patterning geometry, specified, for example, 
in the form of the electrode surface potential, (j)s(x, y) = <|)(x, y, z=0), solve the 
Laplace equation to obtain the potential and the electric field distribution, E (x, y, z) 

30 = -V(t)(x, y, z); for several relevant configurations, analytical models of (t)s(x, y) = 
<|)(x, y, z=0) are available (see Example 1); 

2. Evaluate the term V{E*E); 

3. Evaluate the electrokinetic fluid velocity, of Vf""^"" = -e/uV(|), at 
35 the sUpping plane of the electrode; this serves as the boundary condition for the 
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o evaluation of the flow field above the electrode, further assuming the fluid to be 
static at infinite distance fi-om the electrode; 

4. Determine K(fi)) using the expression: ^(co) =27ia^eM(es(®)- 
eM)/(es(G))+2eM); 

^ 5. Determine the particle velocity, dx^dt = Vb(t) + Vf(t), in the 

particle moving frame; the maximum (relative) particle velocity, Vmax> is attained in 
the vicinity of steep impedance gradients; determine the particle position, Xb = x\i{t); 

\J, J^Ol«^l 11114 IV ^UOIllg UIW lUllVtlVllO VFt/LCUllWU 111 OtUpa -T CUIU ^ UAW 

particle velocity as a function of position: dxb/dr -J{x\i), 
^® In certain embodiments of the present invention, the algorithm detailed in 

the flow chart presented in Fig. 5 is used to convert the experimentally determined 
values for the relaxation frequency, ©r, and the fi-equency-dependent maximum 
particle velocity, Vmax, attained near impedance gradients (see e.g. Fig. 6) into the 
15 particles' surface ("zeta") potential, ^s, and surface conductivity, a. As indicated in 
Figs. 3 and 7, both o>r and Vmax depend on surface potential and ionic mobility and 
the determination of a pair of parameters ((|)s, u) such that cor = cor''^ and 
Vmax = Vnujx is Carried out using numerical analysis as summarized in the flow 
2Q chart. The present invention also includes a method for particle analysis that allows 
simultaneously determination of morphological parameters such as size and shape 
along with the electrical parameters in a miniaturized imaging format. 

This method provides the basis for a quality control technology that is 
particularly sensitive to surface functionalization. In addition, the method provides 
the basis to detect the binding of biomolecules to "matching" molecules displayed 
on the bead surface. That is, the change in surface chemical composition that results 
from the binding of certain molecules to particle or cell surfaces is detected without 
the need for fluorescent or other optical labeling. 
30 The present invention will be better understood from the Experimental 

Details that follow. However, one skilled in the art will readily appreciate that the 
specific methods and results discussed are merely illustrative of the invention as 
described in the claims which follow thereafter. 
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o EXAMPLES 

General Experimental Conditions: 

Preferred particle array assembly is carried out in fluidic microcells which 
are formed by sandwiching a double-sided Kapton spacer of --100 \xm thickness 
5 between a 1 cm x 1 cm substrate that can be patterned electrically, such as a silicon 
chip (n-type, capped either by a uniform or a lithographically patterned thin Si02 
layer), also serving as the bottom electrode, and a glass cover slip coated with 
indium tin oxide (ITO) to a typical sheet resistance of 1400 Ohm Square serving as 
the top electrode. AC voltages of 1-10 Vp.p, in a frequency range from lOO's of hertz 
to a few kilohertz are applied between the electrodes across the fluid gap. Particle 
transport and assembly were monitored by video microscopy permitting frame 
capture and digitization of frames for ftirther analysis. 

Aqueous Suspensions - Prior to assembly, particles stored in buffer are 
15 washed and resuspended in deionized and ultrafiltered water (conductivity < 50 S 
cm*^); in some cases, salt (< ImM) is added to provide buffering. Anionic and 
cationic particles from 0.5 |im to 10 j^m in diameter, composed of silica or modified 
polystyrene and functionalized with a variety of chemical surface groups, as well as 
2Q functionalized core-shell particles obtained from a variety of manufacturers have 
been used. These particles were further modified to display oligonucleotides, DNA 
fragments or proteins. 

Non- Aqueous Suspensions - Illustrative of these systems are: 
DMSO(eDMSO=47.24): 4.4 jam carboxylated PS particles (Dyno, Oslo, Norway); 2.1 

25 

[xm carboxylated PS particles stained with fluorescein (Molecular Probes, Eugene, 
OR); in both these cases, stock particles originally suspended in water/buffer, were 
resuspended in DMSO by way of several successive cycles of centrifugation and 
resuspension to exclude water. Glass beads (3-10 |Lim in diameter) (Polyscience, 
30 PA); stock particles, supplied in dry powder form, were resuspended in the solvent 
of choice prior to assembly, Acetonitrile (eacetonitriie 36.64): 1 jam nonporous oxirane 
acrylate beads (Sigma, MO); procedure analogous as in the case of glass beads. 
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o Example 1- Particle Array Assembly: Half-space Geometry 

Fig. 6 displays bead velocities measured in the course of array assembly of 
oligo-dT(25) particles of diameter 2.8 ^un (obtained from Dynal, Oslo, Norway) 
using Ka = 45 (see also the conditions shown in the inset of the lower panel) and 
^ invoking video microscopy. Analysis of the "streaking" in the snapshots in the 
upper portion of Fig. 6, taken with 1/30 s exposure, reveals the significant 
acceleration of particles traversing the impedance step, created by a variation in 
oxide thickness from ^1000 A to '-1 00 A. Visual inspection also reveals that 
particles are simultaneously displaced toward the substrate. The instantaneous 
10 velocities extracted by image analysis of a series of snapshots are plotted in the 
lower panel of Fig. 6. Arrays of assembled particles act in the manner of low 
mipedance regions of the substrate: they generate in-flow directed toward the array 
perimeter and attract additional particles which accelerate toward the array. 
15 The numerical solution of the equation of motion, given an impedance step 

in a "half-space" geometry, reveals a strong dependence of v^ax on surface potential 
for several ratios of u/uco and for given particle size, as illustrated in Fig. 7. This is 
as expected: a high surface potential corresponds to a large induced dipole moment 
and hence to a large ponderomotive force; in contrast, the dependence on ion 
mobility is weak. 

The equation of motion takes into account the ponderomotive forces, the 
drag forces and also the moving frame of reference of the fluid. The path of any 
bead can be simulated given the experimental conditions such as voltage, gap size, 
25 ionic strength, ion mobilities, oxide thicknesses and given pertinent bead properties 
such as size and surface potential. In general, it is desirable to apply the numerical 
model to fit experimental data and so to improve the accuracy of the determination 

of Vmax. 

2Q A typical solution of the equation of motion, generated numerically by 

representing the "half-space" geometry by a surface potential profile in the form 
(|)s = i4 tanh(icc), is shown in Fig. 9. The origin of the graph represents the location of 
the impedance step. The position of the bead is negative when it is in the thick oxide 
region, where the velocity of the bead is determined by the fluid velocity at any 
given point. However, when approaching the impedance step to within a distance of 
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o the order of the Debye length, the bead experiences the aforementioned 

ponderomotive force which causes it to accelerate across the step relative to the 
fluid. Given the short range of ponderomotive forces (--VJF^), the acceleration ceases 
within a short distance of the step, and the bead experiences a drag force which 

^ decelerates it to the velocity of the surrounding fluid. At larger distances from the 
step, the lateral fluid motion ceases and consequently, the bead comes to rest unless 
the ponderomotive forces surrounding pre-existing bead arrays or clusters provide 
ail additional acceleration to ward Uie array bouiidary. 

Example 2 -Dependence of Relaxation Frequency, cor, on Particle Size and 
10 Composition 

Fig. 10 displays the results of determining the relaxation frequency by 

inducing an array expansion for a variety of particle sizes and for different types of 

particle surface chemistries. In a first group, particles included were carboxylate- 

modified or amine-modified polystryrene particles (obtained from Dyno, Oslo, NO) 

as well as '*oligo-dT(25)" particles (obtained from Dynal, Oslo, NO). In a second 

group, particles were coated with a monoclonal antibody against interleukin-6 

(obtained from R&D Systems, Minneapolis, MN). The anticipated scaling behavior 

of the relaxation frequency, cor^^^ -^a (particle size), as discussed herein, is 

20 

confirmed by the experiments. 

The solid lines represent a fit of the data to a simple expression discussed 
herein and produces a diffusivity parameter, D, which reflects the local ion- 
mediated surface conductance, as discussed herein. Invoking Eq. 4, the data in Fig. 
25 10 may be interpreted to indicate that, for given a and k, the binding of proteins, or 
other large molecules to the bead surface reduces the surface conductivity, a e^nlfy 
introduced herein. 

Example 3 - Particle Analysis 

30 From experiments illustrated in Figs. 6 and 10, values for the maximum 

velocity of beads traversing the step from thick oxide to thin oxide, 
^niax ~ 0.8 mm/s, and the relaxation frequency, cOr =^ 2 kHz, were determined. The 
theoretical model and algorithm underlying Figs. 3, 5 and 7, yield corresponding 

35 values of - 0.5V and u/u^ - 3.9 x 10"^ ; this combination in turn yields a value of 
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a a - 0.15 S/m and Du = 290; here, m„ - 0.02 S /C, as determined from the 
tabulated molar conductivity for Na*. The value for may be compared to the 
surface potential of --330 mV corresponding to the titratable charge of carboxylate- 
modified polystyrene particles of 1 COO 74 or 400 \xC/cm^; the results obtained 
indicate that the effective surface charge of oligonucleotide-modified particles 
substantially exceeds that of carboxylate-modified particles. The value for ©r 
indicates the surface mobility of small cations relevant to the experiments to be 
smaller by a factor of >100 than the mobility in the bulk; however, the surface 
conductivity still exceeds the bulk conductivity by two orders of magnitude, 
reflecting the high value of the surface potential. 

Example 4 - Fractionatioii of Particle Mixture Using Size-Dependence of (Dr 

Fig. 1 1 illustrates the fractionation of a binary mixture of smaller (2.1 ^m 
diameter, characteristic frequency (Ok) and larger (4.5 \im diameter, characteristic 
frequency cor^) particles using a surface design in which a region of lower 
interfacial impedance is created by a variation in oxide thickness from -1000 A to 
-100 A within the stripe indicated by lighter shading (see also Example 1). When 
the frequency, co, of the applied voltage adjusted such that cdr^ < co < (Dr^, only the 
smaller particles experience an acceleration toward the low impedance region where 
they assemble into an array. By narrowing or widening the stripe, flow along the 
stripe also may be generated and may serve as the basis for two-dimensional 
fractionation in the cross-flow so realized. 

Example S - Fractionation of Particle Mixture Using Composition-Dependence 
of OOr 

Fig. 12 illustrates the fractionation of a binary mixture of red and blue 
particles of identical size but differing surface chemical composition. The relaxation 
frequency, cor***"*, of carboxylate-modified particles, tagged with a green dye, is 
lowered by binding of protein, tagged with a red dye, to a characteristic frequency, 
cor'"*^. Under the condition cor"*"'* < co< a)R*'^'^ only the unmodified particles enter 
the stripe region of lower impedance, while modified, protein-coated particles align 
on the high impedance sides of the feature. The magnitude of the shift in 
characteristic frequency is determined by the amount of bound protein. 
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o Example 6 - Detection of Protein Binding to Beads or Cells wittiout the Use of 
Tags 

A molecular binding event on a known particle surface causing a change in 
the surface chemical property (i.e. a change in the dielectric property of the particle 
plus bound molecule over that of the particle) can be easily be detected utilizing the 
5 fact that the effect of the binding is reflected in change in the relaxation frequency 
of the particle. For example, Figure 10 illustrates the experimentally evaluated 
change in relaxation frequency of a particle which display, on its surface, either 
carboxylate groups or carboxylate plus captured layer of antibody. The binding of 
jQ antibody to the particle surface lowered the characteristic frequency, ©mod, of 
modified particles relative to the characteristic frequency, cobarc, of unmodified 
particles. Similarly, the binding of biotinylated oligonucleotide probes on the 
surface of a streptavidin fiinctionalized particle may be detected. The binding of the 
oligonucleotide fragment increased the characteristic frequency, comod, of oligo 
bound particles relative to the characteristic frequency, ©bare, of unbound particles. 



15 



Example 7 - Frequency-Dependent Collection and Assembly of Human 
Myeloid Cells 

Fig. 13 illustrates the collection and enrichment of HL-60 human myeloid 
20 cells (CCL-240, obtained from ATCC and maintained for experiments in 8.5% 

sucrose/0.3% dextrose) in designated areas of a patterned silicon substrate. Using 
silicon chips with patterned oxide surfaces featuring a stripe of thin oxide (< 100 A) 
typically 100 -200 urn in width delimited by ''steps" in impedance indicated by the 
25 change in brightness, application of a voltage of 5 V peak-to-peak, with frequency 
set to the values indicated in the two panels in the figure, produces aggregates of 
either adjacent to, or within the designated stripe, reflecting the selection of the 
frequency in the applied voltage. 

30 Example 8 - Assembly of Bacterial Cells and Fractionation by Selective 
Adhesion 

Fig. 14 illustrates the collection and enrichment of E.coli and 
Staphylococcus epidermis bacterial cells, serving as respective model systems for 
Gram-negative and Gram-positive bacteria, in different regions of a patterned 
35 substrate (here a silicon chip). That is, using chips with patterned oxide surfaces 
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o . featuring a strip of thin oxide (< 1 00 A) typically 1 00 -200 ^m in width, 

apphcation of a voltage of 5 V peak-to-peak across a fluid gap of 100 \im and 
frequency set to 2 kHz produces "bands" of cells in positions reflecting the 
preferential adhesion of S. epidermis near the "step" in impedance indicated by the 

^ change in brightness and E.coli in the center of the strip. The position of the bands is 
thus characteristic of bacterial type. 

Example 9 - Strength of Adhesion 

Cell rolling along a surface is the initial step in many kinds of inflammatory 
response. The present invention provides for the experimental determination of the 
strength of adhesion between cells of interest or ligand-functionalized 
microparticles disposed on a receptor-fanctionalized electrode surface. For example, 
given a patterned electrode surface as described, low-impedance regions can coated 
with E-Selectin to investigate the assembly of HL-60 cells (see example 7). Under 

1 5 the action of LE APS-mediated flow, cells, while being transported toward low- 
impedance regions of the electrode surface, exhibit rolling and sticking, the latter 
representing a manifestation of a positive adhesive interaction between cells and 
surface. In-vitro assays have been described which use the rolling of cells or ligand- 

2Q functionalized microparticles on a receptor-functionalized surface that is exposed to 
tangential fluid flow to estimate the strength of adhesion between cells or 
microparticles and the surface [Brunk, D.K., and Hammer, D.A., Biophysical 
Journal, 72,1997, 2820-2833; Hammer, D.A., and Apte, S.M, Biophysical Journal, 
63,1992, 35-57; Springer, T. , Cell, 76, 301, 1994.]. The analysis of rolling and 
sticking of HL-60 cells on E-Selectin coated glass wafers under the presence of 
pressure driven flow has been used for anti-inflammatory drug screening [Mossavi 
et al. Proceedings of MicroTAS '98 Workshop, D. J. Harrison and A. van den Berg 
(Eds.), 69-72, 1998]. The present invention provides a superior mode of flow 

30 control by way of LEAPS. For example, a receptor-functionalized electrode on 

which cells or microparticles of interest have been disposed in predesignated areas 
can be exposed to LEAPS-mediated fluid flow (see Fig. 4) whose direction and 
strength are readily controlled. This permits the evaluation of the correlation 
between the strength of the flow-induced shear stress acting on the adsorbed cells or 

35 

particles and the efficiency of the hydrodynamic detachment, which in turn permits 
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an estimation of the adhesive forces operative between the cell or particle and said 
surface. 

Example 10 - Multi-dimensional Fractionation 

In multi-dimensional fractionation, a mixture of particles is separated into 
subpopulations in two or more stages in accordance with multiple discriminating 
parameters. Here, in addition to the frequency dispersion of field-induced 
polarization other attributes such as differential selective adhesion etc. can be used. 
For example, a mixture of HL-60, E. coli and S. epidermis can be fractionated using 
patterned chips featuring a strip of thin oxide (< 100 A) typically 100 --200 [im in 
width. Application of a voltage of 5 V peak-to«peak across a fluid gap of 100 nm 
and frequency set to 2 kHz produces "bands" of cells in positions reflecting the 
preferential adhesion of S. epidermis near the "step*' in impedance and E.coli in the 
center of the strip. Under this condition the HL-60 cells form linear chains at the 
thick-thin oxide boundary which extend to the thick oxide region ( see example 7). 

Example 11- Sequential Light-directed Fractionation of Multi-component 
Mixtures 

A multi-component particle mixture in which each particle population 
displays a unique relaxation frequency (cor ) may be fractionated using a 
combination of a patterned electrode and illumination as illustrated in Figure 15. 
Given relaxation frequencies (cor ) in the range ©low < o)r < ©high, the frequency of 
the applied AC electric field is scanned from the high limit (cohigh) to the low limit 
(cDiow) of the range; particles are sequentially swept into the thin oxide/low 
impedance collection region starting with the ones with higher relaxation 
frequencies. Once inside, each population is guided by light into separate thin oxide 
pads which are connected to the initial collection zone via narrow channels. This 
particular patterning will confine beads in the respective pads. 

Example 12 - Encoding of Particle Arrays using Engineered Particle 
Polarizability 

Dielectric "tagging" to produce particles of desired polarizability and 
corresponding relaxation frequency may be used to produce a specific dependence 
of the particle response to LEAPS. Dielectric tagging may be accomplished by way 
of introducing labels to the exterior or interior of the particle. For example, 
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o metallization of particles changes their polarizability. This approach may be 

extended to the engineering of a particle whose polarizability is light dependent. For 
example, using light-sensitive tags such as Ti02 nanoparticles, the polarization of a 
Ti02-tagged particle is increased exclusively when illuminated with light [Y. 

^ Komoda, Tata N. Rao, and A. Fujishima, Photoelectrorheology of Ti02 nanoparticle 
suspensions. Langmuir, 1997, 13, 1371-1373]. 

Although the invention has been described above with reference to examples 
and to nreferred embodiments, it will be anoreciated that the invention mav be 
embodied in other forms or carried out in other ways without departing from the 
10 spirit or essential characteristics thereof. The above description is therefore to be 

.considered in all respects illustrative and not restrictive, and all changes which come 
within the meaning and range of equivalency are intended to be embraced therein. 
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